Microscopic analysis of the valence band and impurity band theories of (Ga,Mn)As 



O 



r-- 



J. Masek, 1 F. Maca, 1 J. Kudrnovsky, 1 0. Makarovsky, 2 L. Eaves, 2 R. P. Campion, 2 K. W. Edmonds, 2 

A. W. Rushforth, 2 C. T. Foxon, 2 B. L. Gallagher, 2 V. Novak, 3 Jairo Sinova, 4 ' 3 and T. Jungwirth 3 ' 2 

institute of Physics ASCR, v.v.i., Na Slovance 2, 182 21 Praha 8, Czech Republic 

2 School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, United Kingdom 

3 Institute of Physics ASCR, v.v.i., Cukrovarnickd 10, 162 53 Praha 6, Czech Republic 

4 Department of Physics, Texas A&M University, College Station, TX 77843-4242, USA 

We analyze microscopically the valence and impurity band models of ferromagnetic (Ga,Mn)As. 
We find that the tight-binding Anderson approach with conventional parameterization and the full 
potential LDA+U calculations give a very similar picture of states near the Fermi energy which 
reside in an exchange- split sp-d hybridized valence band with dominant orbital character of the host 
semiconductor; this microscopic spectral character is consistent with the physical premise of the kp 
kinetic-exchange model. On the other hand, the various models with a band structure comprising an 
impurity band detached from the valence band assume mutually incompatible microscopic spectral 
character. By adapting the tight-binding Anderson calculations individually to each of the impurity 
band pictures in the single Mn impurity limit and then by exploring the entire doping range we find 
that a detached impurity band does not persist in any of these models in ferromagnetic (Ga,Mn) As. 

PACS numbers: 75.10.Lp,75.30.Hx,75.50.Pp 
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Over more than four decades, (Ga,Mn)As has evolved 
from a pioneering direct gap p-doped semiconductor [1 
into an archetypical degenerate semiconductor with hole- 
mediated ferromagnetism [2-4 . Paramagnetic insulating 
Gai-^Mn^As materials prepared in the 1970's by melt 
growth showed valence band (VB) to impurity band (IB) 
activation, with a non-systematic filamentary metallic 
conduction being observed at the highest studied dopings 
of x ~ 0.1% and ascribed to sample inhomogeneities [5]. 
The degenerate semiconductor regime was not reached in 
these materials. A comprehensive experimental assess- 
ment of basic doping only trends became possible in the 
late 1990's with the development of epitaxial (Ga,Mn)As 
films [2j [7]-[9] which can be doped well beyond the equi- 
librium Mn solubility limit while avoiding phase segrega- 
tion and maintaining a high degree of uniformity. Trans- 
port measurements on such films confirmed the insulating 
characteristics and the presence of the IB for x < 0.1%. 
For higher concentrations, 0.5 < x < 1.5%, no clear sig- 
natures of activation from the VB to the IB have been 
detected in the dc transport, suggesting that the bands 
start to overlap and mix, yet the materials remain insu- 
lating. At x ~ 1.5%, the low-temperature conductivity 
of the films increase abruptly by several orders of magni- 
tude and the material becomes a bulk degenerate semi- 
conductor. The onset of ferromagnetism occurs on the 
insulating side of the transition at x ~ 1% and the Curie 
temperature gradually increases with increasing doping, 
reaching ~190 K at the accessible substitutional MnQ a 
doping of x ~ 8%. 

One physical scenario for ferromagnetic (Ga,Mn)As, 
termed the VB picture, has an exchange-split band struc- 
ture comprising the impurity band merged into the va- 
lence band. The states at the Fermi energy, Ep, retain 
the predominant orbital character of the host semicon- 



ductor and are moderately hybridized with the localized 
Mn d-electrons. This description, quantified by a va- 
riety of theoretical methods, has been a fruitful basis 
for analyzing and predicting a whole range of thermo- 
dynamic, magnetic, transport, and optical properties of 
ferromagnetic (Ga,Mn)As [3j[4]. Recently, several exper- 
imental observations have been interpreted using alterna- 
tive models of an IB which is detached from the VB [TOl - 
H4] . However, it has been argued that a detailed analysis 
of the data in combination with transport experiments is 
also consistent with the VB picture [7 . The postulated 
IB models have not been previously derived from a mi- 
croscopic theory considering all relevant orbital states in 
the mixed crystal. In order to help resolve the debate on 
these alternative interpretations, we examine here the IB 
models by recreating them using microscopic modeling 
techniques and studying their band structure character- 
istics over the entire doping range. These calculations (i) 
firmly establish the microscopic basis and internal con- 
sistency of the VB picture, (ii) demonstrate the mutual 
inconsistency of the various postulated IB models, and 
(iii) demonstrate that a detached IB does not persist 
in any of the models' band structures at dopings corre- 
sponding to ferromagnetic (Ga,Mn)As. Our theoretical 
analysis is based on the tight-binding Anderson (TBA) 
approach which includes all spectral components in the 
band structure forming the states near Ep, accounts for 
the Mn ^-orbital electron-electron interaction effects us- 
ing the self-consistent unrestricted Hartree-Fock method, 
and can be adopted to realize microscopically the diverse 
proposed IB models. Additional physical insight is pro- 
vided by comparisons to full-potential LDA and LDA+U 
calculations. More details on the techniques and more 
extensive numerical results can be found in the Supple- 
mentary material [T5] . 



The perturbation of the crystal potential of GaAs due 
to a single Mn impurity has three components, (i) The 
first is the long-range hydrogenic-like potential of a sin- 
gle acceptor in GaAs which produces a bound state at 
about 30 meV above the VB |B]. (ii) The second contri- 
bution is a short-range central-cell potential. It is spe- 
cific to a given impurity and reflects the difference in the 
electro-negativity of the impurity and the host atom [16] . 
For a conventional non-magnetic acceptor Znca? which is 
the 1st nearest neighbor of Ga in the periodic table, the 
atomic p- levels are shifted by ~0.25 eV which increases 
the binding energy by ~ 5 meV. For Mn, the 6th near- 
est neighbor of Ga, the p- level shift is ~1.5 eV which 
when compared to Znc a implies the central-cell contri- 
bution to the acceptor level of MnQ a ~ 30 meV [17] . (iii) 
The remaining part of the MnQa binding energy is due 
to the spin-dependent hybridization of Mn d-states with 
neighboring As p-states. Its contribution, which has been 
directly inferred from spectroscopic measurements of un- 
coupled MnQa impurities [I7J US] , is again comparable 
to the binding energy of the hydrogenic single-acceptor 
potential. Combining (i)-(iii) accounts for the experi- 
mental binding energy of the MnQ a acceptor of 0.1 eV. 
An important caveat to these elementary considerations, 
further quantified by our microscopic calculations [15] , is 
that the short-range potentials alone of strengths inferred 
in (ii) and (iii) would not produce a bound-state above 
the top of the VB but only a broad region of scattering 
states inside the VB. 

The VB picture of ferromagnetic (Ga,Mn) As builds on 
the above conventional semiconductor description of the 
MnQa acceptor in which the presence of the long-range 
hydrogenic-like impurity potential is essential for creat- 
ing a bound state in the band gap. With increasing dop- 
ing, the impurity level broadens and for a sufficiently 
screened hydrogenic potential the impurity states must 
merge into the VB within this picture. The premise of the 
models with a persistent detached IB in ferromagnetic 
(Ga,Mn)As is distinct and can be reconciled by ascribing 
the main role in binding to the short-range potentials and 
a minor role to screening and impurity level broadening. 

We now provide microscopic analysis of these scenarios 
by performing TBA band-structure calculations. Disor- 
der is treated in the coherent potential approximation 
(CPA) which allows us to scan the entire range of dop- 
ings from the single Mn-impurity limit to MnAs. Our re- 
sults are consistent with available corresponding spectra 
obtained using the super-cell method [19] which justifies 
the validity of the CPA [15] to represent the one-particle, 
orbital resolved, density of states (DOS). We first take 
the conventional parameterization of atomic levels and 
overlap integrals [15j [16]. On-site electron-electron in- 
teractions on the Mn d-states are described using the 
Hubbard parameter U — 3.5 eV and the Heisenberg pa- 
rameter Jh = 0.6 eV, which also correspond to a conven- 
tional parametrization of d-orbital correlations in atomic 



Mn or Mn in II- VI semiconductors, and are consistent 
with values of U and Jh inferred from photoemission 
experiments in (Ga,Mn)As [20]. Since we are primarily 
interested in the ferromagnetic behavior which occurs at 
relatively high dopings (> 1%) and is governed by the 
spin-dependent p-d hybridization potential, we omit the 
long-range Coulomb potential which is non-magnetic and 
largely screened at the relevant hole concentrations [21]. 
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FIG. 1: (a),(b) Total (solid line) and partial Mn d-orbital (red 
filled) DOSs calculated using the TBA (CPA) and compared 
with the LDA+U (super-cell) results for Mnc a dopings of 
3,6,12, and 100%. (c),(d) Same for the TBA d compared with 
the LDA results. Ef is indicated by the dashed vertical line. 

Using the conventional values of the TBA parameters 
[15J[T6] we fi rs t determined occupation numbers on Mn d- 
orbitals and corresponding on-site energies using the self- 
consistent unrestricted Hartree-Fock description of the 
Anderson Mn impurity embedded in the semiconductor 
environment [15 . The important result of these calcu- 
lations is that we did not find any tendency to symme- 
try breaking in these occupation numbers, i.e., the three 
^2^-orbitals (and similarly the two weakly hybridized e g - 
orbitals) remain degenerate and strongly localized. After 
determining the Mn <i-orbital on-site energies we proceed 
to calculate the microscopic DOS of (Ga,Mn)As over the 
entire doping range. In Figs. l(a),(b) we plot examples 
of both the total (black line) and the Mn d-orbital re- 
solved (red filled) DOSs for x =3, 6, and 12% together 
with the results for x =100%, i.e. for the zinc-blende 
MnAs. The Mn d spectral weight is peaked at approxi- 
mately 4 eV bellow the top of the VB, in agreement with 
photoemission data [20] , and is significantly smaller near 
Ef as further highlighted in Fig. 2(b). The Fermi level 
states at the top of the VB have a dominant As(Ga) p- 
orbital character; the stronger As p-component is plotted 



in Figs. 2(a). The p-d coupling strength, 7V /3 = A/(Sx) 
[3], determined from the calculated VB exchange split- 
ting A (and taking S = 5/2) is close to the upper bound 
of the reported experimental range of Nq/3 ~ 1 — 3 eV 
[T71I2Q1I22H 21, as shown in Fig. 2(c). This is regarded as 
a moderately weak p-d coupling because the correspond- 
ing Fermi level states of the (Ga,Mn)As have a similar 
orbital character as the states in the host GaAs VB. The 
spectral features shown in Figs. l(a),(b) and 2(a)-(c) are 
among the key characteristics of the VB picture. Note 
that the kp kinetic-exchange (Zener) model calculations 
assume a value of 7V /3 also within the range of 1-3 eV 
(typically closer to the lower experimental bound) |3]. 
It is this moderate p-d hybridization that allows it to 
be treated perturbatively and to perform the Schrieffer- 
Wolff transformation to effective valence band states ex- 
periencing a spin-dependent kinetic-exchange field [3]. 
Hence, the effective kinetic-exchange model and the mi- 
croscopic TBA theory provide a consistent physical pic- 
ture of ferromagnetic (Ga,Mn)As. 
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FIG. 2: Partial As p-orbital (a) and Mn d-orbital (b) DOSs 
at Ef in the depicted TBA models. The remaining contri- 
bution to the total DOS is primarily due to Ga and Mn p- 
orbitals with their relative weights given nearly precisely by 
the Ga/Mn ratio [15]. (c) The p-d coupling strength in the 
TBA models compared with the LDA (triangles) and LDA+U 
(circles) results, (d)-(f) Total DOSs showing the merging of 
the IB into the VB in the depicted tight-binding models. 

We next attempt to recreate the IB models by consid- 
ering that the bound state at the single Mn is formed 
by the short range impurity potentials. As noted above, 
this is not obtained from the conventionally parameter- 
ized TBA and the values of the atomic levels or overlap 
integrals have to be adjusted ad hoc to match the 0.1 eV 
binding energy [15] . We first search for a bound state 
due to the central-cell potential by treating the Mn p- 
level as a free parameter. We find that binding the hole 
by the central-cell potential alone requires physically in- 



comprehensible p-level shifts of several tens of eV [25] . 
The reason for this is the short-range nature of the po- 
tential and the orbital composition of the top VB from 
which the bound state forms. The VB near its maximum 
is dominated by p-orbitals of the As not Ga sublattice. 

A more favorable scenario to create a bound state 
through the short-range potentials is tuning the strength 
of the p-d hybridization. This term is less local as it af- 
fects four As neighbors of the Mnc a impurity and acts 
on the As p-orbitals which form the top of the host VB. 
To tune the hybridization strength we can treat as a 
free parameter the atomic Mn d- level or the p-d hopping 
[T5] , The corresponding models are labelled as TBA d 
and as TBA pd , respectively. For the TBA d model we ob- 
tain the 0.1 eV bound-state when shifting the d-level by 
1.5 eV. We now fix this parameter and calculate the cor- 
responding DOSs over the entire doping range, as shown 
in Figs. 1(c), (d). More detailed characteristics of the cor- 
responding spectra are summarized in Figs. 2(a)-(d). The 
key observation is that for dopings above ~ 0.1% the 
band structure cannot be recast in a model with Fermi 
level states residing in a narrow IB (of width not ex- 
ceeding the single impurity binding energy) which is de- 
tached from the VB. The TBA pd model yields the same 
general conclusion, as shown in Figs. 2(d),(e). A de- 
tached IB model for the ferromagnetic (Ga,Mn)As ma- 
terials is therefore microscopically incompatible with the 
0.1 eV acceptor level even if the binding of the hole to 
the Mn impurity was entirely due to short-range poten- 
tials. We remark that Nq/3 in the TBA d ( pd ) parameteri- 
zation is about a factor of 2-3 stronger than in the con- 
ventionally parameterized TBA model, i.e., much larger 
than the upper experimental bound for the p-d coupling 
strength. This discrepancy is due to the omission of the 
long-range Coulomb potential when fitting the experi- 
mental single MnQa acceptor state. Note also that the 
dip in the TBA d ( pd ) DOS, which persists and shifts deep 
in the band at high dopings, is another consequence of 
the ad hoc increased p-d coupling. 

We next associate the IB models postulated in lit- 
erature with corresponding microscopic TBA calcula- 
tions. One proposed phenomenology assumes a domi- 
nant Mn ^-orbital nature of the detached IB and allows 
for some hybridization with the host VB [10] [14] . The 
TB A d theory is the closest microscopic realization of this 
model. It shows that apart from the absence of the de- 
tached IB itself at dopings corresponding to ferromag- 
netic (Ga,Mn)As, the formation of the 0.1 eV acceptor 
state by shifting the Mn d-level does not yield a dominant 
Mn d spectral weight near Ep. From this perspective we 
regard the 0.1 eV acceptor level as moderately shallow. 

An orthogonal IB model, in terms of the assumed or- 
bital character of the IB, elaborates on a sp-tight-binding 
Hamiltonian with shifted ^-levels on four As neighbors 
of the Mn Ga [13 . We label this model as TB*\ The 
shifts are introduced to effectively account for the micro- 



scopic p-d hybridization and again to obtain the 0.1 eV 
single impurity state without the hydrogenic long-range 
Coulomb potential. The model has a merit in the very 
dilute regime [25 as the extent of the bound state wave- 
function (the exponential tail) is determined by the value 
of the binding energy and is insensitive to the specific 
choice of the confining potential. It also captures, by 
its design, the symmetries of the As p-orbital dominated 
bound state. The model can be associated with our mi- 
croscopic TBA pd calculations and indeed the correspond- 
ing DOSs show very similar doping trends, as shown in 
Fig. 2 (e),(f). Again, no detached IB persists in the TB P 
DOS to dopings above ~ 0.1%. 

Another phenomenological proposal assumes that 
states in the IB have Mn p-orbital character [11]. This 
corresponds to our first attempt to obtain the MnQ a ac- 
ceptor level in short-range potentials by considering the 
central-cell component only. As discussed above, such 
a model would require an unphysical large shift of the 
Mn p-levels. If we omit the microscopic justification of 
this IB model the approach has a merit as a phenomeno- 
logical effective model describing the Mn acceptor level 
in the band-gap of the host semiconductor. Since the 
impurity states are added in this type of effective model- 
ing ad hoc to the spectrum, the model does not conserve 
the total number of states. (In a microscopic language 
it describes rather an interstitial than a substitutional 
impurity.) The applicability of the approach is therefore 
limited to small Mn concentrations and the model is not 
suitable for exploring trends with changing Mn doping. 

Finally we show in Figs. l(a)-(d) examples of the com- 
parisons of the TBA d and TBA calculations with results 
of the LDA and LDA+U full-potential ab initio theo- 
ries [15j 125] . We find a very good agreement between 
the TBA and LDA+U results [27]. The LDA+U, the 
TBA, and the kinetic-exchange Zener theories therefore 
all provide a compatible picture of the band structure 
of ferromagnetic (Ga,Mn)As. We also find a clear cor- 
respondence between the TBA d and LDA results. The 
large exchange splitting of the VB obtained in the LDA 
reflects the general deficiency of the LDA to account for 
localized states within an itinerant band. Mn ^-states 
are more delocalized and move closer to the VB edge in 
the LDA, which enhances the hybridization. Hence, the 
exchange splitting of the top of the VB is increased to 
values comparable to those of the TBA d Hamiltonian. 

To conclude, at the doping levels for which (Ga,Mn)As 
is ferromagnetic none of the postulated one-particle DOS 
models with a detached IB arising from the 0.1 eV Mn 
acceptor level in GaAs is microscopically justified. The 
Fermi level states in ferromagnetic (Ga,Mn)As can be 
regarded as residing in a modified VB of the host semi- 
conductor due to disorder, exchange splitting, and ad- 
mixture of the impurity orbitals. The corresponding 
one-particle band structure can be described by meth- 
ods ranging from full-potential density-functional the- 



ory to multi-orbital tight-binding- Anderson or envelope- 
function approaches which are all mutually consistent. 
We also emphasize nevertheless, that due to the vicinity 
of the metal-insulator transition and correlation phenom- 
ena these effective one-particle VB band models can only 
represent a proxy to the complex electronic structure of 
ferromagnetic (Ga,Mn)As materials. 
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SUPPLEMENTARY MATERIAL 



TIGHT-BINDING ANDERSON DESCRIPTION OF (Ga,Mn)As 

In this section we review our implementation of the tight-binding Anderson formalism for calculating the band 
structure of (Ga,Mn)As. First we introduce the Hamiltonian describing the semiconductor host GaAs. We use the 
Slater-Koster tight-binding approach [I] [2] which is a variant of the linear-combination-of-atomic-orbitals (LCAO) 
method implemented using the Bloch's theorem. Most of the calculations in this paper are done in the second- 
nearest-neighbor (nearest Ga-As, As- As, and Ga-Ga bonds) sp 3 parametrization [3] of the tight-binding Hamiltonian. 
We have checked that the alternative nearest-neighbor (only nearest Ga-As bonds) sp 3 parametrization [4] and the 
nearest-neighbor sp 3 s* parametrization [5 give very similar results for the physics discussed in this paper. Note that 
these parametrizations are obtained by fitting to the band structure calculations with pseudopotentials obtained from 
experiment. Anticipating the introduction of substitutional Mnc a impurities we formally included Ga d orbit als in 
the GaAs tight-binding Hamiltonian. There energy level is, however, far above the bottom of the conduction band 
and has no effect on the part of the band structure we are interested in. For the description of a single Mn impurity we 
use a combination of the Bloch's theorem based Slater-Koster approach in which Mn does not introduce new orbitals, 
i.e., it does not change the dimension of the Hilbert space, electron-electron interactions on the Mn d shell are treated 
in the Anderson-impurity spirit, and the effect of the periodic crystal environment of the host GaAs on Mn impurity 
states is accounted using Green's function formalism. 

Slater-Koster tight-binding description of GaAs in the second-nearest-neighbor sp 3 parametrization 

For describing the valence and conduction bands of the host GaAs, the basis functions are written in the form of 
Bloch sums, 

N-l / N-l \ 

*kaa = JV-VV*-*' J2 eik ' R ^aa(r -R n -p a )=l N -l/2 e ik-p ae ik.r ^ g-^^^r -Rn- Pa ) 
n=0 \ n=0 / 

N-l 

= N' 1 ' 2 J2 e zk ' Pna 0aa(r - p na ) with p na = R n +p a , (1) 

n=0 

where a is the atom index in the unit cell, a is the atomic orbital quantum number, N is the number of atoms a (unit 
cells), n is the unit cell index, R n is the unit cell vector, p a is the position vector of the atom a in the unit cell, and 
the overall phase factor e lk ' Pa is added to obtain more symmetric expressions for the Hamiltonian matrix elements 
derived below. The wavefunction (j) aOL (r — Pna) are centered around the atom a in unit cell n and are orthonormal, 
i.e., assumed to be constructed from the atomic orbitals i/j aa following Lowdin's orthonormalization procedure, 

(t>aa{r ~ Pna) = ^aa(r - p na ) - - ^^ aQ (r - pma)Smn + • • • , (2) 

m 

where S mn are the overlap integrals, 

dr ^la( r ~ Prna)lpaa(r ~ Pna) (3) 



/- 



For the zinc-blende GaAs lattice we set, 



p a = (0, 0, 0) for Ga 

= ^(1,1,1) for As (4) 

where a\ c is the cube lattice constant of GaAs. 

The number of k points in the first Brillouin zone is the same is the number of unit cells but since the Hamiltonian 
has the same periodicity as the basis function in Eq. ([I]) it is diagonal in the k- vector. 

(§k>a>cx>\H\<& kaa ) = — ^ ^ / rfr<ft*/ a /(r ~ Pn'a')H(/) aa (r - p na ) 

n n' ^ 

n J n' 

J2e ikPna J dr^ aa (r - Pna )H^^r) (5) 



N 



In this Koster-Slater implementation of the tight-binding model the problem is simplified as we only diagonalize a 
Hamiltonian matrix of a dimension given by the number of atoms in the unit cell and valence orbitals considered for 
each atom. Since every atom in GaAs has the same environment (there's one Ga and one As per unit cell), As-Ga 
Hamiltonian matrix elements, keeping only the nearest-neighbor (4 atoms) hopping terms in the sum, can be written 
for all k- vectors within the first Brillouin zone as, 



3 n.n. ^ 

H Asa , Gal3 (k) = J2 eik ' PnA ' / dr<t>* Asa (r - p„ A8 )ff<W) • (6) 



The unit cell lattice vectors are a\ = ^(1, 1,0), a 2 = ^(1,0, 1), and a 3 = ^(0, 1, 1). For the Ga in the 0th unit 
cell considered above, Rq = (0,0,0), the 4 nearest As neighbors are in unit cells i^o, Ri = — ai, R2 = —Q>2, and 
R3 = — as at positions, 

POAs = P As + #0 = -^(1,1,1) 

PlAs = PAs + Rl = -^(-1,-1,1) 

p2As = PAs + #2 = -y (-1,1,-1) 

PSAs = PAs + R3 = ^(1, "I, -1) (7) 

The Ga-Ga matrix elements, keeping the diagonal term and only the nearest Ga (12 atoms) hopping terms in the 
sum, read 



12 n.n. „ 

H Gaa , Gal3 (k) = J2 eikR " / dr ^haa(r - Rn)H<p GaP (r) , (8) 

n=0 ^ 



with the positions of the 12 As neighbors given by, 



Hi = -ai = ^(-1,-1,0) 

R2 = -a 2 = ^(-1,0,-1) 

R3 = -a 3 = ^(0,-1,-1) 

^4 = ai = ^(1,1,0) 

#5 = a 2 = ^(1,0,1) 

He = a 3 = ^(0,1,1) 

R 7 = a 2 -03 = -y (1,-1,0) 

# 8 = a 3 -a 2 = — -(-1,1,0) 

# 9 = ai -a 3 = — -(1,0, -1) 

#io = «3 -01 = — -(-1,0, 1) 

#11 = ai -a 2 = -y (0,1,-1) 

R 12 = a 2 -ai = — -(0, -1,1) 



Similarly the As- As elements read, 



The diagonal in the atom index on site integrals are denoted as, 



£ a = J dr<j> Gaa (r)H<j> Glia (r) 
t s = ydr0Ua(»WABa(r) 



The diagonal in the atom index hopping integrals are parametrized as, 



The off-diagonal in the atom index hopping integrals are parametrized by 

^ S ' Ga (PnAs) = / dr(f)* Asa (r - pnAs)H(l) Ga(3 (r) 



(9) 



12 n.n. „ 

H AsaMP (k) = J2 eikRn / dr <t>l s «( r - Rn)Hc/> ABP (r) . (10) 

n=l ^ 



(11) 



WgCfln) = J dr<p Gaa (r - R n )H<f> Ga p(r) 

W^(R n ) = J dr4>* Asa (r - R n )H4> AsP (r) . (12) 



(13) 



Finally we rewrite the hopping integrals in the form of a and tt bonds between s, p, and d orbit als as, 



W^(R n ) - V^QRnl) 



W£ a pj (Rn) 

W s A s S (Rn) 

W£ S pj (Rn) 

W^ Ga ( PnAs ) 

w£> G *( PnAa ) 

W^ Ga ( PnAs ) 

w A p f a ( PnAs ) 

W^f*( PnAs ) 
W^ y2 ( PnAs ) 
W^ 2 a _J PnAs ) 

<S a (P«As) 
W p % a (PnAs) 

W Px d 2 2 (PnAs) 



W PyC\^PnAs) 



W-As,Ga / \ 

PzdoAPnAs) 



W-As,Ga / \ 

P,d 3z2 _ r2 (PnAs) 

W-As,Ga / \ 

P,d 3z2 _ r2 (PnAs) 



V G *(\R n \)R nti ,i = x,y,z 



V^RnDR^iRnj + V^{\R n \)(8ij - R n ,iR n ,j) ,i,j =x,y,z 

V s t(\Rn\) 



= V s %(\Rn\)Rn,i ,i = x,y,z 



Vj^dRnDR^iRnj + V^dRnDiSij - R n>i R nJ ) ,i,j = x,y, 

V s A s ^(\ PnAs \) 
Kpa Ga (|PnA s |)/5„A s> i ,i = x,y,z 

^a Ga (|PnAs|)/3„AM ,i = X,y,Z 
ppcr \ 

V^Vsda a (\PnAs\)PnAs,iPnAs^ 



V pp* a (\PnAs\)pnAs,iPnAs,j + V pp % a (|p n As|)(^j ~ PnAs,iPnAs,j) ,hj = X,y, 

As Cla .,, u. . . . 

,hJ =x,y,z 



V ^ T/ As,Ga/i i\ / -2 -2 \ 

~^~ V sda \\ Pn As\) [PnAs,x ~ PnAs,y) 



V s A d* Ga HPn As |) 



^2 



PnAs,z cy 



\PnAs,x ' PnAs,y) 

^ V pda *(\PnAs\)P n As,iPnAs,j + V^ Ga (|p nA s|)PnAs,j (l ~ 2p£ AM ) , 

vSVpda a (\PnAs\) Pn As,iPn As J PnAs,k ~ 2V p dn ^PnAs \) PnAs,iPn As, j PnAs,k ,^j^ 

V^ T rAs,Ga/i i\a / -2 -2 \ , TrAs,Ga/i i\- /-, -2 , -2 

~2~ V pda \\PnAs\)PnAs,x [PnAs,x ~ PnAs,y) "•" V^ UPnAs|;PnAs,x (^ ~ PnAs,x "•" PnAs,^ 

~2~Vpda \\PnAs\)PnAs,y [PnAs,x ~ PnAs,y) ~ Vpdn (\ Pn As\ )Pn As, y (1 + P n As,x ~ PnAs,y 

V^TrAs^a/i i\ - / -2 -2 \ TrAs,Ga/i i\ - / -2 -2 \ 

~^~ V^r \\PnAs\)PnAs,z {PnAs,x ~ PnAs,y) ~ V pdir {\PnAs\)PnAs,z [PnAs,x ~ PnAs,y) 



V pd* a (|PnAs|)/3nAs,^ 



V pda a (\pnAs\)pnAs,z 



-2 
PnA 



s.z 2 

1 



(,PnAs,a; i PnAs,y) 



-2 _ / -2 | -2 \ 

PnAs,z ~ 9 \PnAs,^ ' PnAs,y) 



~ ^ V p A dn Ga (\PnAs\)pnAs,ip 2 nA s,z , i 
+ V^^'^dPnAsDPnAs^ (PnAs,* + 



a, 2/ 



-2 \ 

PnAs,y) 



(14) 



Here p n As,i (-Rn,i)is the z-th component of the unit vector along p n As (Rn)- Note that the diagonal in the atom 
index hopping integrals are parametrized within the three-center approximation, i.e., considering also the term in the 
periodic potential, V(r) = ^Z na V a (r — R n — p a ) centered between the wavefunctions on site corresponding to the 

other sublattice. In this approximation we need to keep spatial index (x,y,z) dependent parameters V Pi p i7T (\R n \). 
For the near-neighbor As-Ga hopping we only consider the two-center integrals in which case all V's including the 
^?p?r' Ga (|PnAs|) term depend only on the bond length. Values for the on-site energies and all the above a and 7r-bond 
hopping energies for s and p orbitals are taken from Ref.|3l For the Ga d orbitals, which are unimportant for GaAs but 
are included anticipating the Mn substitution, we put the on-site energy to a very high value. The hopping energies 

T/As.Ga T rAs,Mn T rAs,Ga T rAs,Mn i T rAs,Ga T rAs,Mn j. j • j_i r n • -l- 

V sda = V sda > V pda = V pda and V P dn = V pd* are discussed in the following section. 

We have tested that spin-orbit coupling, when included in the diagonal Ga and As p-orbital terms in the form of 
^Ga(As)^ • s ? nas no significant effect on the physics discussed in this paper. (Note that in the sp 3 s* parametrization 
[5], spin-orbit is included in the fitting of the tight-binding parameters to the ab initio band structure. The second- 
nearest-neighbor sp 3 parametrization [3] was done in the original paper without including spin-orbit coupling; in our 
implementation, ^Ga(As) were taken from Ref. [5] and the on-site energies were then rigidly shifted to recover the value 
of the GaAs band gap.) 



Tight-binding parametrization of a single substitutional Mnc a impurtity in GaAs 

As mentioned above, the GaAs tight-binding Hamiltonian parameters [3H5] are obtained by fitting to the empirical 
pseudopotential band structure calculations. For the nearest-neighbor s and p-orbital As-Mn hopping energies we 
take the same values as for the As-Ga bonds. The Mn <i-orbitals is included in our theory in the spirit of Harrison's 
rules [2 j. The hopping integrals involving Mn <i-orbitals are parametrized as, 

T/ As,Mn _ ^ frd\ 3/2 

T/ As,Mn _ & frd\ 3/2 

v p d* ~ r hd^ med 2 \ d ) 

T/ As,Mn & (r d \ 3 / 2 ,,-, 

y Pd , =^^UJ ' (15) 

where we took for the dimensionless hopping parameters for the covalent Mn-As bonds, rj s da — —4.9964, rjpdo- — 
—4.6644, and r\ v d^ = 2.1503, the Mn-As distance d = \/3/4 a^ c , and for the effective Mn d-orbital radius Td = 0.86A 
(see Ref. [2j). The change of the crystal potential due to the substitutional Mn is represented by replacing Ga on- 
site energies with on-site Mn energies, e^ n , e^ n , and e^ s at the site where the substitution takes place. This is 
a natural tight-binding representation of the central cell correction [7]. For the on-site Mn d energies, the t<i g 
(m = xy, xz, yz) and the e g (m = x 2 — y 2 , 3z 2 — r 2 ) orbitals are split by the crystal field, A c f. The other reason why 
we keep explicitly the d-orbital index m and the spin index s = ± is that the levels can be further spit by electron- 
electron interaction. This effect can be modeled by the multi-orbital Anderson many-body Hamiltonian solved in the 
unrestricted Hartree-Fock approximation. The direct and exchange Coulomb terms are parametrized by 



U mm , = J dr J dr'r ms (rW m 'Ar')V(\r-r'\)^ ms (r)4> m , s ,(r') 

Jmm' = [dr f dr'r ms (rWm' S (r')V(\r - r'\)^ ms (r')4> m , s (r) . (16) 



In the unrestricted Hartree-Fock theory the occupation numbers, (n m>s ), of the d-orbitals hybridized with the host 
GaAs are not a priori know and have to be calculated self-consistently This procedure will be explained below; here 
we just assume for the moment that they are known. Using the above parametrization, the Hartree-Fock equations, 

$Z,.il>m.(r) = (E d + A m )x[, ms (r)+Y, f dr '\ip m 's'(r')\ 2 V(\r-r'\)xp ms (r) 

m' ,s' 

- J2 I ' dr'r m >s(r')ip ms (r')V(\r - r'|)Vw s (r) , (17) 



yield 



<%Z,s =Ed + ^m+Yl U ™m> ((nm' t s') ~ f) ~ J2 J ™™' H 71 ™' ,*) ~ f) ' ( 18 ) 



Here Ed is the d- level of an isolated atomic Mn calculated from ab initio assuming unpolarized Mn atom, i.e., 
(nm,s) — f — 1/2, A m = 2/5 A c f for the three t^g orbitals and A m = — 3/5 A c j for the two e g orbitals, / = 1/2 is 
subtracted from the occupation numbers in the Hatree-Fock terms to recover Ed in the unpolarized configuration of 
an isolated Mn atom. 

Anderson [8] reduced the number of direct and exchange interaction parameters by setting U mm ' = U for all m, m' 
and J mm' = J for all m ^ m' '. The interaction Hamiltonian in this parametrization is, however, not rotationally 
invariant in the spin space which can be corrected [9j [10] by considering U mm ^ U mm r, namely U mm t = U for 



m ^ m! and U mm = (U + J). (Note that /7 mm = J mrn (see Eq. (16) which removes the self-interaction.) We use this 
parametrization in which the on-site Mn d energies read, 

e™% = E d + A m + (U + J) «n m ,_ s ) -/) + ([/- J) J2 « n «^> - /) + U E ««m'.-) " /) • ( 19 ) 

m'y^m m'y^m 

The values of the relevant parameters are summarized in Table I. 
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As 
c s 
Ga 

Mn 


-6.724 
-3.978 
-0.200 


As 

Ga 
t p 

Mn 


0.641 

2.874 
4.874 


E d 
Ac/ 
U 
J 


-0.100 
0.500 
3.500 
0.600 


V s da 
Vpda 
ypdix 


6.770 

-6.320 

2.913 


Mn 
6 d t 2g,+ 

Mn 
€ d eq ,+ 


-2.220 
-3.013 


Mn 
6 d t 2g,- 
Mn 

d e g ,— 


2.634 
2.364 



TABLE I: Parameter of our tight-binding model for Mn in GaAs. All energies are given in electronvolts and referred to the 
top of GaAs valence band parametrized in Ref. [3] 

The self-consistent occupation numbers are determined by employing the local Green's function formalism. The 
formalism allows us to evaluate local and orbital resolved density of states on the Mn in the environment of the 
GaAs host. We remark that the same formalism can be used to evaluate the Mn d-orbital occupation numbers in 
the environment of (Ga,Mn)As with the Green's functions describing the (Ga,Mn)As environment obtained from the 
coherent potential approximation (discussed below). The important result of these occupation number calculations 
is that they remain very similar in the GaAs and (Ga,Mn)As environment and in both cases not very far from the 
polarized isolated Mn configuration. More details on this are given later in these supplementary notes section; in the 
following section we recap the formalism of the calculations. 



Green's functions formalism for evaluating Mn d-orbital occupation numbers 



As mentioned above, <i-orbitals are also included in the tight-binding description of the cation sites occupied by the 
host Ga atoms. Atomic levels of Ga d-orbitals are, nevertheless, chosen far from the valence band so that their mixing 
with the band states is practically excluded. In this way we can use the established sp 3 second-nearest-neighbor 
tight-binding parameters of GaAs, while having a tight-binding basis consisting of 26 orbitals per unit cell (18 spd 
orbitals on the cation site and 8 sp orbitals on the anion site). In this basis the introduction of the Mn impurity 
does not require to enlarge the Hilbert space but is only represented by shifted on-site energies of the spd orbitals 
(Koster-Slater impurity [6 ); recall that the hopping energies are considered to by the same for Ga and Mn. 

The Hamiltonian of the GaAs host crystal in the representation of the basis Bloch sums (Eq. (fil)) can be written as 



H (k)=D + W(k), 



(20) 



where D° is a diagonal matrix (26x26) with Ga and As on-site energies on its diagonal, and the "kinetic energy" 
matrix W(k) is constructed from the hopping integrals Ga-As, Ga-As, and As- As hopping integrals. 
The propagator (retarded Green's function) from one atom to another in the GaAs lattice is given by, 



k>a>(Pn>a')\Go(t)\(j) aa (p na )) = (<j> a > - a > < (p n > V )|e %hH °\(j>aa{pna)) 



In Fourier transform, 



G (z) 



1 f°° 

— / dte* z G(t) = (zt-H )- 1 , 
in Jo 



where z = e + irj. In the representation of the basis Bloch sums (Eq. (HI)) we obtain, 

G (z) = Yl \®ka'a>)($kaa\ (zl - JfoW)^,, 

k(B.z.),a'a / ,aa 



(21) 



(22) 



(23) 



The orbital-diagonal on-site elements of the Green's function for the 0-th unit cell on the Ga sublattice then read, 

G , a ( Z ) = (4>Gaa(0)\Go(zMGaa(0)) = ^ £ (zl - H^k))^ = fi u . c . J j^ (zt - ffflW)^ , (24) 



k(B.z.) 



(2tt) 2 
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where fl u . c . is the unit cell volume and the fe-sums and integral are over the 1st Brilloun zone. We neglect the 
orbital-off-diagonal on-site elements of the Green's function. (Note that due to the tetrahedral symmetry, most of the 
orbital-off-diagonal on-site elements of G(z) are zero.) 

Neglecting the long-range Coulomb part of the impurity potential of the substitutional MnQ a , we can describe Mn 
by a Koster-Slater model |B]. The Koster-Slater model assumes that the perturbation is restricted to the impurity 
site and represented by shifts of the atomic levels. For the MnQa impurity we write 

8 a ,s = %!-£*. (25) 

The Hamiltonian of GaAs with a single Mn impurity at site "0" then reads, 

H(k) = H {k)+H 1 , (26) 

where 

Hi = y^5 a , s \a)(a\ 

\a) = \<pGaa(0)) (27) 

The Green's function for H can be written as, 

G = (zl-H - iJi)" 1 = ((zl - H )(l - (zl - Ho)- 1 ^))- 1 = G (l - Go^i)" 1 

= Go J2( G oH!) j =G + GoH.Go + G Q H X G Q H X G Q + ... (28) 

3 

Because of the orbital-diagonal form of both the on-site Green's function matrices and the perturbation, the calculation 
of the impurity Green's function reduces to a set of scalar equations. We obtain, 

G ajS = G 0j a + ^2(a\Go\a f )S a f i3 (a f \Go\a) + ■ ■ ■ 

Oi' ,s 

— Go^ a + Go,q;Oq;,sGo,q; + G() ya a )S Gq )(x a )S G Q )(x . . . 

— Go^ a + Go,q;Oq;,sGo,q; / (^a^Co^J 

3 

Go ' a (29) 



1 — 3a,sGo, a 

Finally we rewrite the Green's functions in a physically and computationally convenient way, 

GoAz) = (z-£a a -Fa(z)y 1 , (30) 



and (29 (yield 



which, together with Eq, (24), should be regarded as a definition of T a (z). This definition useful because Eqs. p0 



G atS (z)=(z-e™ n s -r a (z)y 1 , (31) 

where the "self-energy" T a (z), calculated from the host GaAs Green's function, represents the effect of hopping 
energies and is the same in G as in Go, i.e., does not depend on the shifted energies on the impurity site. In other 



words, Eq. (31) shows how the local electronic structure of the impurity depends both on the impurity potential (ej^) 
and on the electronic structure of the surrounding crystal (T a (z)). 

From the orbital-diagonal on-site Green's functions we obtain orbital decomposition of the local density of states 
(DOS), 

ffa,,(e) = --ImG a , a (e + iO) . (32) 

7T 

The desired expression for occupation numbers of the Mn d-orbital states are then given by, 

(nm,a) = / degm iS (e)f EF (e) , (33) 
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where the Fermi function /# F (e) = 0(Ep — e). The self-consistent unrestricted Hatree-Fock calculation of the Mn 
d-orbital on-site energies, e^ n s , start from assuming (n m ^) = 1 — 5 m?+ for the five majority spin d-orbitals and 
{rim-) = (5 m> _ for the five minority spin d-orbitals. The initial small and d-orbital dependent variations from 1 
(or 0) are set to allow for spontaneous symmetry breaking which could result, e.g., in only four Mn d-electrons 
remaining strongly localized deep in the band and the fifth d electron state shifted towards or above the top of the 
host semiconductor valence band. The self-consistency loop then proceeds by putting these initial occupation numbers 



in Eq. (19) for the Mn on-site energies, then calculate the impurity Green's function from Eq. (I3l| and obtain new 
occupation numbers from Eqs. ( [32] ) and ( [33] ). These self-consistent calculations can be performed assuming pure GaAs 
environment, i.e., using the function T a (z). The function T a (z) can be also replaced with T^ s (z) which describes 
the environment of (Ga,Mn)As in the coherent-potential approximation (CPA). The CPA method is described in the 
following subsection. 

Coherent-potential approximation 

At finite concentration of Mn the single-impurity picture breaks down as soon as the typical distance of the 
impurities becomes comparable with the extend of the related wave functions. We evaluate the band structure of a 
mixed (Ga,Mn)As crystal using the (CPA), which is particularly suitable for describing the system over the entire 
doping range from to 100% of Mn. The basic quantity in the CPA is the configurationally averaged Green's function. 
The averaging restores the translational invariance of the mixed crystal. The procedure is based on replacing real 
but random atomic levels by site-independent but complex and energy-dependent selfenergies S aQ; , s which are to be 
determined selfconsistently. Because of the diagonal form of the Green's functions, S aa)S will be also diagonal and 
since there's in disorder in (Ga,Mn)As on the As-sublattice, 

X Asa = et s ■ (34) 

The self-energy on the Ga-sublattice £c7aa,s is obtained by from the CPA condition which states that the on-site 
Green's function for a periodic system in which the Ga-sublattice energies have been replaced by ^Gaa,s = ^a,s equals 
the composition- weighted sum of the single-impurity on-site Green's function assuming Ga or Mn energies on the site: 

G 0,as( Z ) = C 1 ~ %)GG^a,s( Z ) + ^Mn^,^) • (35) 



Here x is the doping in the Gai-^Mn^As mixed crystal. Using Eqs. (24) and ([30]) we can write the Green's function 
of the effective periodic Hamiltonian, with the diagonal terms Dq in Eq. (|20|) replaced with £ acM , as 



Gl a Jz) = {z-E a , s (z)-Tl s (z)) 1 ■ (36) 

Similarly the on-site impurity Green's functions read, 

G Ga , a , s (z) = (z-e^-Tl^z))' 1 

G M n,aAz) = (s-^-lt.C*))" 1 . (37) 



The convergence of the self-consistent solution of Eqs. ( |36| ) and (37) is guaranteed by starting the calculations from the 
the virtual-crystal approximation which assumes that E a?s = T> v a c s are the average atomic levels on the Ga sublattice 
for the given doping, 

S- = (l-x)6^ + <^. (38) 

RESULTS: (Ga,Mn)As BAND-STRUCTURE CALCULATED USING THE UNMODIFIED TBA 

PARAMETRIZATION 

Cartoon representation of the expected band structure of Mn doped GaAs 

Fig. HI a) shows the expected ~ 30 meV binding energy due to the long-range Coulomb potential of a single acceptor 
in GaAs. Fig. [lib) shows the comparable enhancement of the binding energy due to the central cell correction whose 
strength is estimated form the conventionally parametrized on-site energies on Mn and from comparison of the 
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a) 



long-range Coulomb 



-30 meV 




c) 



-0.1 eV Mn Ga acceptor state 



short-range p-d hybridization 



short-range central cell 

m 





b) 



short-range central cell 



-30 meV 




Mn Ga " acceptor 



d) 



short-range p-d hybridization 



no bound-state above V.B. 




broad scatt. states in V. 
Mndt 




FIG. 1: Schematic diagrams of the long-range (a) and the two short-range (b),(c) contributions to the single Mn acceptor level. 
(d) No bound state in the gap would form without the long-range Coulomb potential. 

Mn(Ga) substitution with the well established central cell correction for the non-magnetic Zn(Mn) substitution (see 
main text and references therein). Fig.fijc) shows the again comparable enhancement of the binding energy due to the 
hybridization of the Mn d-states with the host valence band states (see main text for references). Fig. [TJd) illustrates 
that no bound state above the top of the valence band is expected to form in the absence of the long-range Coulomb 
potential of the charged Mn acceptor. When this potential is screened in the high doped metallic ferromagnetic 
(Ga,Mn)As, it is then plausible to expect that also no detached impurity band can persist above the valence band 
and that the two must merge. 



Single Mn impurity: Local density of states and occupation of Mn d-orbitals 



The self-consistent unrestricted Hartree-Fock occupation numbers calculated for a Mn impurity in a pure GaAs 
environment (T a (z)) can be classified according to their spin and to the tetrahedral symmetry of the impurity, 
i.e., there are three degenerate ^-orbitals and two weakly hybridized e^-orbitals. The electron interactions within 
Mn d-shell are strong enough to result into almost complete spin polarization ((n t2 +) = 0.899, (n eg ^+) = 0.979, 
(n t2g -) = 0.108 and (n e _) = 0.038) and to formation of a local magnetic moment close to 5/i#. On the other hand, 
we did not find any tendency to the quenching of the orbital moment, i.e. symmetry breaking in these occupation 
numbers due to the electron correlation. [8 This is a key observation allowing us to disregard a formation of Mn-related 
bound state due to the correlation effects and stick to an effective one-particle picture of the impurity state. 

Fig. El shows spin-polarized local density of states (LDOS) at a single substitutional Mn impurity in GaAs in 
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FIG. 2: Local density of states at Mnoa and its orbital decomposition. 



the spectral range of valence and the lowest conduction bands. The shaded area in the left panel represents the 
contribution of Mn d-states. The change of local densities of s and p-states at the Mn impurity, together with the 
corresponding local densities of states at the host atom Ga, are shown in the right panel. Their differences reflect the 
central-cell corrections to the atomic levels e s and £ p , respectively. 

An alternative picture of the impurity induced modification of the electron spectrum is obtained from the change 
Sg tot (E) of the total DOS in which the changes of the local densities of states at all lattice sites are summed up. The 
integrated quantity 



SN tat (E) 



f 

J — c 



Sg tot (E')dE'. 



(39) 



is particularly suitable to show spectral features with a small weight at the impurity site. In our case of Mno a , 5N tot (E) 
combines two features, i.e., addition of the d-orbitals and reconstruction of the band states due to both central-cell 
corrections and hybridization. Fig. [3] shows smooth, but non-monotonic increase of SN to t(E) in the valence band. 
This is a signature of mixing of Mn d-orbitals with the band states. In addition to the changes in the valence band, 
a sharp dip in SN tot (E) for both spin polarizations around 2 eV indicates that Mn d-orbitals hybridize significantly 
also with GaAs conduction band states (away from the r-point).[TTJ[T2] 

On the other hand, as the Fermi energy remains pinned at the edge of the valence band in the single impurity 
regime, 5N to t(^F) m the band gap define the number of electrons accumulated at and around Mn impurity. These 
numbers are five and zero for spin-up and spin-down electrons, respectively, as a result of the fact that the impurity 
potential does not pull any state from the valence band in the conventionally parametrized TBA with self-consistent 
Hatree-Fock Mn d-orbital on-site energies when the long-range acceptor Coulomb potential is not included. 
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FIG. 3: Integrated change 5N to t(E) of the total density of states induced by substitution of single Mn for Ga in GaAs crystal. 
Electron states with spin parallel (up) and antiparallel (down) to the local magnetic moment are treated separately. 

DOSs of (Ga,Mn)As over the entire concentration range 

We now fix the Hartree-Fock d-levels obtained for a single Mn impurity and perform the CPA calculations of the 
DOS in Mn doped GaAs as well as in (Ga,Mn)As mixed crystals with higher concentrations of Mn. We show the 
atom and orbital resolved DOSs for the entire doping range up to MnAs. Fig. [4] shows the result for GaAs doped 
with 10 percent of Mn. The admixture of Mn <i-orbitals to the unoccupied states at the top of the valence band is 
relatively small and the occupation numbers of the <i-orbitals do not differ significantly from the values obtained in 
the single impurity case ((^t 2ff ,+) = 0-886, (n e +) = 0.984, (nt 2 -) = 0.110 and (n e ) = 0.035). This confirms that 
also the self-consistency condition for Mn atomic levels remains fulfilled with a reasonable accuracy. 

For higher concentrations of Mn (see Figs. [5][8]), the spectral weight near the Fermi energy gradually evolves from 
being dominated by the As p-orbitals to a comparable As p and Mn d-orbital weight for Mn content around 70 
percent (see Fig. [7|) where the DOS undergoes a transition from the shape characteristic to Mn-doped GaAs to the 
form similar to the MnAs limit (see Fig. |8|. 
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FIG. 4: TBA density of states of Gao.gMno.i As and its orbital composition. Position of Fermi energy is indicated by a vertical 
line. 
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FIG. 5: TBA density of states of Gao.7Mno.3As and its orbital composition. Position of Fermi energy is indicated. 
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FIG. 6: TBA density of states of Gao.5Mno.5As and its orbital composition. Position of Fermi energy is indicated 
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FIG. 7: TBA density of states of Gao.3Mno.7As and its orbital composition. Position of Fermi energy is indicated 
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RESULTS: (Ga,Mn)As BAND-STRUCTURE CALCULATED USING MODIFIED TBA 

PARAMETRIZATIONS 

Cartoon representation of the modified TBA parametrizations 



a) 

0.1 eV 



^Sb 



c) 



Mn-d 



^Z^b 



short-range p-d hybridization 




short-range p-d hybridization 




Mn eft Aspt 



b) 

0.1 eV 



Mn-p | 



^Sb 



d) 



As-p* 



^ZS5 



short-range central cell 




short-range p-d hybridization 

m 



Mndt 




Mndt Aspt 



FIG. 9: Cartoon representation of the modified TBA parametrizations. 

Fig. KM a) illustrates the strategy of calculations in this section in which we modify the conventional TBA parametriza- 
tion in order to obtain the 0.1 eV bound state of a single Mn impurity in GaAs without the long-range Coulomb poten- 
tial. Fig.[9jb) shows the first, central-cell approach in which we shift the Mn p-orbital on-site energies. Figs. [9jc),(d) 
show the enhanced p — d hybridization approach in which we shift the Mn <i-orbital on-site energies or increase the 
values of the hopping energies involving Mn d-orbitals. 



Enhanced central-cell correction by shifted Mn p-orbital on-site energies 



Here we consider modified on-site Mn p-orbital energies, 



-Me 



Mr 



a Mn 



(40) 



We examine a capability of such perturbation to create a localized state in the band gap by using the impurity Green's 
function, Eq. [3l] A bound state is represented by a pole of the Green's function, i.e., the energy z of the bound state 
is obtained from a condition 



-Mn 



r P (*) 



(41) 
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FIG. 10: Relation between the p-orbital energy e^ n at substitutional Mn impurity in GaAs and position of related bound 
state in the band gap. Observed Mn-related acceptor level is represented by horizontal dotted line. The arrow indicates the 
conventionally parametrized TBA value of e^ n . 

Fig. [lO] shows that no bound state can be formed in the band gap of GaAs for any physically meaningful value of 

a Mn 



Enhanced p — d hybridization by increased Mn d-orbital on-site energies — TBA d 



Here we consider modified on-site Mn d-orbital energies obtained by shifting the one-particle energies in Eq. (19) 



as, 



E d 



k Mn 



(42) 



Fig. 11 shows that A^ n « 1.6 eV is necessary to create a bound state with a binding energy 0.1 eV above the edge 
of the valence band. Fig. 12 shows the transformation of the impurity related resonance in the valence band into a 
bound state in the band gap for increasing value of A^ n , indicating the same value of A^ n which gives the 0.1 eV 
bound state. 

This value of A^ n is not as unphysical as in the above case of shifted Mn p-levels. However, it represents a transfer 
of Mn (i-orbital energies from « -2.5 eV to ~ -1 eV and, as a result, a significant enhancement of the hybridization 
effects leading, e.g., to a factor of 2 enhancement of the exchange coupling J p d between the local moments and the 
holes. This value of J p d is much larger than any experimentally inferred value of this parameter. 

We now use the modified TBAd with A^ n = 1.59 eV that results in the localized state at 0.1 eV in the single 
impurity case and use the CPA to visualize the evolution of the impurity band with increasing concentrations x of Mn. 



Fig. 13 shows that the width of the impurity band quickly increases with increasing x. The impurity band merges 
with the valence band already for x < 0.1% and at x ~ 6% the width of the broad maximum of the deformed valence 
band extends over the lower half of the band gap. 



Finally, the left panel of Fig. 14 shows the spectral density for the mixed system with 4% of Mn in the spectral 
range close to the top of the valence band. We note that around the center of the Brillouin zone, the spectral density 
has a double-maximum structure. The lower maximum corresponds to the host valence band, the upper one to the 
scattering states at the top of the valence band due to Mn impurities. This maximum broadens and finally disappears 
for wave-vectors away from the center the Brillouin zone. The effective mass corresponding to the scattering and host 



parts of the band are similar as indicated in the right panel of Fig. 14 . 
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FIG. 12: Integrated change SNtot(E) of the total density of states for spin-up electrons visualizes formation of a bound state 
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FIG. 13: Total density of states of dilute (Ga,Mn)As for A^ n = 1.59 eV adjusted to the position of the impurity level to the 
observed acceptor binding energy 0.1 eV. 
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FIG. 14: Spectral density of states of Gao.96Mno.04As for Ajf n = 1.59 eV along L-r-X path in the Brillouin zone. The maxima 
of the spectral density are used to show effective dispersion relations in the right panel. 
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Enhanced p — d hybridization by increased Mn d-orbital hopping energies — TBA pd 



Here we consider modified Mn d-orbital hopping energies obtained by multiplying the r] coefficients in Eq. ( 15 ) by 
a factor A > 1, 



Vpda — -^]pd<j 
Vpdir = -A-VpdTt - 



(43) 



Fig. [15| shows that A ra 1.5 eV is necessary to create a bound state with a binding energy 0.1 eV above the edge 
of the valence band. Similarly to the TBA P case, J p d is enhanced by a factor of 2 in TBA pd . To show what can be 
expected at finite concentrations of Mn, we adopted the modified hopping integrals with A = 1.5 and calculated the 
density of states shown in Fig. 16 Even though the mechanism of formation of the bound state was different, its 



transformation into a broad spectral feature merged with the valence band is very similar to the results from Fig. 13 
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FIG. 15: Relation between the energy of Mn- induced related bound state in the band gap of GaAs and the enhancement factor 
A for sp-d hybridization. Dotted line represents Mn acceptor level and the arrow corresponds to our basic parametrization 
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FIG. 16: Total density of states of dilute (Ga,Mn)As for the sp 
acceptor level at 0.1 eV. 



d hybridization factor A adjusted to obtain the observed 



RESULTS: (Ga,Mn)As BAND-STRUCTURE CALCULATED USING THE sp TIGHT-BINDING MODEL 

WITH SHIFTED As p-LEVELS - TB P 



Here we show results obtained using the conventional tight-binding parametrization of GaAs and an effective model 
of the Mn impurity based on shifting majority-spin As p-orbital on-site energies on the neighbors of the presumed Mn 
impurity. 



-As 

-As 



As 



.As 



k As 



(44) 



First we find and fix the value of A^ s by considering a single As impurity with e^ s and searching for the 0.1 eV 

shows that A^ s « 6 eV must be used. Then we do 
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bound state without the long-range Coulomb potential. Fig. 
the CPA to calculate the band structure as a function of the presumed Mn doping. Before writing down the CPA 
equations we need to realize that the single As impurity problem with shifted As p-orbital energies on a single As 
atom corresponds to 3/4 of one Mn impurity. The factor 1/4 is because Mn interacts with 4 As nearest-neighbors 
and the factor of 3 is because Mn is not distributed symmetrically around the As but the bond is directional. From 
the As p x , p y , and p z orbitals only one of the three independent linear combinations is directed along the bond with 
Mn; for example, for the pQAs = ^(0,0,0) bond, it is the combination p x + p y +p z . The CPA equation then reads, 



G 0,a,s ( z ) = (1 - x)G A s,aA Z ) + xG As,a,s( Z ) 



(45) 



where the single-site Green's function GAs,a,s( z ) corresponds to the unshifted As p-orbitals on the site, e^ s , and 
G As,a,s( z ) t° the shifted orbitals e^, and the effective Mn concentration x — | x where x is the concentration of As 
atoms with shifted p-orbital energies considered in the CPA calculations. 

Also in this case, the width of the broad maximum corresponding to the 



The results are summarized in Fig 18 



Mn-induced scattering states depends on Mn concentration in the same way as for the models studied in the previous 
Section. 
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FIG. 17: Relation between the modified As p-orbital energy e Ps +, representing the exchange field due to Mn impurity in GaAs, 
and a position of related bound state. Dotted line represents Mn acceptor level and the arrow corresponds to the unperturbed 
atomic level e^ s . 
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FIG. 18: Total density of states of dilute (Ga,Mn)As for the A^ s adjusted to obtain the observed acceptor level at 0.1 eV. 



FULL-POTENTIAL SUPERCELL AB INITIO CALCULATIONS IN LDA AND LDA+U 



To support the tight-binding studies, we present the results of the supercell calculations performed by using the full- 
potential linearized- augmented-plane- wave method (WIEN package [13 ). Besides the standard LDA calculations (with 
an improved local-density approximation by the generalized-gradient approximation, the GGA) we used also LDA+U 
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(or GGA+U to be precise) approximation with the same parameters U and J as used in the tight-binding calculations. 
Because the LDA+U approximation is known to shift the occupied ^-states to lower energy and unoccupied ^-states 
to higher energies, the relation between LDA+U and LDA seems analogous to the relation between conventional and 
modified versions of the tight-binding method. 

We present the densities of states calculated for the superlattices GayMnAsg, GaisMnAsie, and Ga3iMnAs32, 
representing the diluted systems with roughly 12%, 6%, and 3% Mn, respectively, and also for a hypothetical MnAs 



crystal with a zinc-blende structure as a limiting case. The results are summarized in Fig. 19 and Fig. pO 
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FIG. 19: FP LAPW densities of states of (Ga,Mn)As supercells representing three typical concentrations of Mn together with 
the DOS of zinc-blende MnAs crystal obtained with LDA (GGA to by precise) approximation. DOSs are normalized to the 
volume of conventional unit cell. 
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FIG. 20: FP LAPW densities of states of (Ga,Mn)As supercells representing three typical concentrations of Mn together with 
the DOS of zinc-blende MnAs crystal obtained with LDA+U (GGA+U to be precise) approximation. DOSs are normalized to 
the volume of conventional unit cell. 



TB-LMTO CPA LDA AND LDA+U CALCULATIONS 



The supercell calculations are restricted to a limited range of chemical compositions. In particular, it is difficult 
to represent very diluted system with Mn concentrations below one percent. We took an advantage of the available 
tight-binding linear muffin-tin approximation (TB-LMTO) version of the CPA [14-16 and performed a series of cal- 
culations for a complementary set of Mn concentrations. Again, we are interested in differences induced by accounting 
the correlation effects and compare the DOSs obtained with LDA and with LDA+U. The results for very diluted 



(Ga,Mn)As magnetic semiconductors are summarized in Fig. 21, densities of states for concentrated mixed crystals 



are shown in Fig. 22 and Fig. 23 
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FIG. 21: Partial density of Mn d-states in dilute (Ga,Mn)As magnetic semiconductors. 
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FIG. 22: Density of states for (Ga,Mn)As mixed crystals with various content of Mn obtained in the LDA approximation. 
Shaded area shows the partial density of Mn d-states. 
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FIG. 23: Density of states for (Ga,Mn)As mixed crystals with various content of Mn obtained in the LDA+U approximation. 
Shaded area shows the partial density of Mn d-states. 
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